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ABSTRACT 

This thesis describes an alternating renewal process and 
a computer program that simulates a transient alternating 
renewal process. The purpose of the simulation is to test 
analytical approximations for the efficiency of a system 
described by such a process. 
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I. 



DESCRIPTION OF AN ALTERNATING RENEWAL PROCESS 



This thesis describes an alternating renewal process and 
gives a computer program that simulates a transient alter- 
nating renewal process. Consider a system that can be in 
one of two states, up or down. The system starts up and 
alternates between up and down, remaining in each state a 
random time. The successive times in each state are jointly 
independent. The sequence of up times are identically 
distributed with distribution and the sequence of down 

times are identically distributed with distribution 

th 

^d* ^i ^i corresponding i^ duration of time 

the system is in the up state and down state. The means of 
the two distributions are assumed to be finite and are 
denoted by and y^ respectively. The up rate and down 
rate are denoted by and X^, and are related to the 
means as 



y = = mean time up 

^u X ^ 

u 



y , = = mean time down 

d X , 
d 



( 1 ) 




Figure 1. An Alternating Renewal Process. 
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Two quantities of interest concerning an alternating 
renewal process are the availability and the efficiency. 

The availability, A(t), is the probability that the system 
will be up at time t. The efficiency, E(t), is the ratio 
of the expected amount of up time during [0,t] to t. The 
relationship between E{t) and A(t), which will be derived 
later, is 

t 

f A(x) dx 

E(t) = (2) 

t 

An alternating renewal process has two imbedded renev/al 

processes. Reference 1 gives a detailed description of a 

renewal process. One process, ca].l it a type 1 process, has 

as its renewals the events when the system goes from the up 

to the down state. This is a general renewal process with 

the first interevent distribution and the remaining 

interevent distributions all being F = F^ o F^, where 'o' 

denotes convultion. The duration of the i^ interevent time 

is for i ^ 2. A second imbedded renewal process 

has as its renewals the events when the system goes up. The 

common interevent distribution is F = F o F,. This will 

u d 

be denoted as a type 2 process. 

The expected number of renewals during [0,t] is denoted 
by M(t) and is called the renewal function. If F is a 
function of bounded variation, that is, the difference of 
two monotone functions, and vanishes to the left of zero, 
then 
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Figure 2. Type 1 and Type 2 Renewals 
and the Interevent Times. 



ic 

F (s) 



/ exp (-sx) dF (x) 
0 - 



is defined as the Laplace-Stiel jes Transform of F. 

times between renewals are distributed as Fj^,F,..., 

* F*(s) 

^ "" 1 - F*(s) * 



If the 
then 



(3) 



For a derivation of this expression see Ref. 1. 

For the type 1 imbedded renewal process with renewal 
function Mj^(t), 



M^(s) = 



F„(s) 



1 - F*(s) F*(s) 



For the type 2 imbedded renewal process with renewal 
function (t) , 

M*(s) = L . 

1 - F*(s) F*(s) 

Define the random variable S(t) such that 

f 1 



S(t) 



2 



if the system is up at t 
if the system is down at t. 



H 





« 




Imbedded Renewal Processes. 



Let = number of renewals of type i, i = 1 , 2 , that occur 
in [0 , t] . Then 

S (t) = - nJ + 1, 



If the expectation is taken, then 



E[S(t)] = E[N^] - E[nJ] + 1 



= M^Ct) - M^(t) + 1 = A(t) 



(4) 



The Laplace-Stiel jes Transform of (4) gives 



A*(s) = M*(s) - M*(s) + 1 



F*(s) F*(s) 



F* (s) 
u 



1 - F*(s) F*(s) 1 - F*(s) F*(s) 



+ 1 



A"(s) = 



1 - F (s) 
u 



1 - F*(s) F*(s) 



(5) 



Let the random variable Q(t) be the amount of up time 
during [0,t] and R(t) be the expected amount of up time over 
the interval [0,t]. Then 
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t 



Q(t) = / S(x) dx . 

0 

If the expectation is taken, then 



R(t) = E[Q(t) ] = E[/ S(x)dx] 

0 

t 

= / A(x) dx 
0 



from which (2) immediately follows 



If (5) can be inverted, then availability and, hence 
efficiency, can be obtained. For example, if the distribu- 
tions of the up and down times are exponential with cor- 
responding rates and A^, then 



t 



E(t) 



t 

/ A(x) dx 
0 



= R(t)/t . 



t 



F^(t) = 1 - exp(-A^t) 



F^(t) = 1 - exp(-A^t) 



( 6 ) 




A . 
1 



u,d 



A . + s 
1 



1 



gives 



A (s) 




+ s 
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Inverting yields 



A(t) = 



Xd+A 






Use of (2) gives 



E(t) = 



/ 

0 



X ,+X 

d u 



A If 



exp[- (X^+X^)x] 



) 



dx 



u 



^u'^^d t(X +X^) 
u d 



2 (1 - exp [- (X^+X^) t] ) 



A (s) is not always easy to invert. For instance, if 
the down distribution is k-erlang, then 



1 - 



u 



a"(s) = 



X +s 

u 



X \ / X, \k 

u I / d i 



1 - 









(A^+s) (Ad+s) - A^A, 



(7) 



The inversion of (7) can be obtained in principle but it is 
not easily expressed in closed form. If the down distribu- 
tion is log normal, a reasonable assumption for many 
systems, then A(t) cannot be obtained in closed form. 

If the system has been operating for a long period of 
time it will settle down into a steady state. For large t 
an asymptotic approximation for E(t) can be found using 
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c? 



and 



lim 

t, 00 



t 

/ A(x) dx 
0 

t 



nl 



lim 
s 0+ 



s 



n-1, * 
A 



(s) 



d F* (s) 
d s 



s = 0 



y 



F*' (s) 



, F is any distribution, 
s = 0 



Letting n = 1, then 



/ A(x) dx 

lim 0 

t +. 



I lim 

II s 0+ 



1 - F (s) 
u 



1 - F*(s) F*(s)J 



Use of L' Hospital's rule for limits in the indeterminate 
form gives 



f A(x) dx 

lim 0 _ lim 

t 4 - s 0 + 



-F*'(s) 

u 



-[F*' (s)F* (s) + F*(s)F*' (s) ] 



u 



y +y j 

'^u d 



( 8 ) 



The derivation of (8) was independent of the type of distri- 
butions. If the system is in a transient state, and equi- 
librium has not been reached, then the asymptotic 
approximation is not valid and other approximations must be 
made. Hence, a need exists for testing the approximations. 
The program described in the next section simulates an 
alternating renewal process for this purpose. 
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II. 



A SIMULATION OF AN ALTERNATING RENEWAL PROCESS 



This section describes a simulation of an alternating 
renewal process. Three distributions are used. 
Exponential 

F(t) = 1 - exp(-Xt) t ^ 0 

= 0 otherwise. 

k-erlang 

k-1 (Bt)^ exp(-Bt) 

F(t) = 1 - I t > 0 

i=0 1* 

= 0 otherwise 



Log normal 



F(t) = / ^ 

0 ax/27T 



exp[- — ~ (S,n x-a)^] dx 
20 ^ 



= 0 



t ^ 0 
otherwise. 



Two combinations of the distributions are used. In one 
combination the up times are distributed exponentially and 
the down times are distributed k-erlang. In the other sit- 
uation, the up distribution is the same but the down 
distribution is log normal. 

R(t) , the expected amount of up time, and E(t), the 
efficiency of the system over the interval [0,t], are the 
output parameters and are provided for ten time intervals. 
Analytical approximations for E(t) are given for comparison. 
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The input parameters vary with the choice of down 
distribution. There are some inputs common to both situa- 
tions: At the unit time interval, the up time rate, 

X^ the down time rate, the number of simulations and runs 
desired. If the k-erlang distribution is selected the 
parameter k must be provided. If the log normal distribu- 
tion is used the variance must be given. The computer 
program gives more information on the inputs. 

The program functions generally as outlined in the flow- 
chart. Each simulation has these basic steps. First the 
up cycle is entered. While in the up cycle and upon 
exiting, checks are made to determine if a time period has 
been spanned. If not, the next step is the down cycle. 

While in the down cycle and upon exiting, tests are made 
again to determine whether or not the end of a time period 
has been reached. If not, the program returns to the up 
cycle. If at any check the time interval over which the 
simulation is to occur has been spanned, and, if the desired 
number of simulations has not been completed, then another 
simulation begins. After the required number of simulations 
have been completed, the output parameters are computed and 
printed . 

The program was written in FORTRAN IV and used on an 
IBM 360/67 computer. 
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E# (t) = the first analytical approximation of E(t) and is the exact data if both 
distributions are exponential. 
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